Inelastic collapse of a randomly forced particle 
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We consider a randomly forced particle moving in a finite region, which rebounds inelastically 
with coefficient of restitution r on collision with the boundaries. We show that there is a transition 
at a critical value of r, r c = e~ 1T ^ i ', above which the dynamics is ergodic but beneath which the 
particle undergoes inelastic collapse, coming to rest after an infinite number of collisions in a finite 
time. The value of r c is argued to be independent of the size of the region or the presence of a 
viscous damping term in the equation of motion. 
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Randomly accelerated particles have been used as 
models in many contexts, ranging from the classic prob- 
lem of Brownian motion ffl to ecological phenomena such 
as swarming H. Theoretical studies have centered on 
first-passage problems both with [Rj and without 
viscous damping or in the presence of a potential J(| . Re- 
cently, systems of such particles that collide inelastically 
have been studied in connection with granular media . 
These systems are characterized by clustering M , which 
is thought to be a collective effect related to inelastic 
collapse of ballistically moving particles Q . 

In this Letter we show that inelastic collapse can occur 
in a system consisting of a single particle. Specifically, we 
study a randomly forced inelastic particle in a finite box. 
If the coefficient of restitution governing the collisions 
with the boundaries is less than some critical value, r c , 
the particle will collide with one wall an infinite number 
of times in a finite time and come to rest at the wall. The 
value of r c is shown to be independent of the system size 
or the presence of viscous damping. This behaviour rep- 
resents a novel example of a broken ergodicity transition 
in a single particle, one-dimensional stochastic system. 
It also has repercussions on the use of randomly-forced 
inelastic particles as models of granular media. 

Consider a particle moving with position x at time t 
within a region of size I and subject to a random force. 
The equation of motion is 
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Hp 
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for < x < I, where rj{t) is a Gaussian white noise with 
correlator (rj(ti)rj(t2)) = 28(ti — t<z) and initially we ig- 
nore viscous damping. When the particle collides with 
one of the walls it rebounds inelastically with coefficient 
of restitution r, i.e., Vf = — ru,;, where Vi and Vf are 
respectively the velocities just before and just after the 
collision. The random force tends to increase the par- 
ticle's energy, while the collisions dissipate energy. One 
might naively expect that the motion would settle into a 
steady state with a well-defined average energy, except in 
the pathological cases r = (where the particle simply 
adheres to the first wall it hits) and r = 1 (where there is 



no dissipation and the particle's energy increases without 
limit). 

Solutions to the Focker-Planck equation for a randomly 
accelerated particle are technically formidable, even in 
time-independent situations where the boundary condi- 
tions are known p[. In the present problem, the bound- 
ary conditions must be determined self-consistently, so 
that the flux of particles leaving the system with veloc- 
ity — \v\ is balanced by a flux of particles re-entering the 
system with velocity r\v\. We have not been able to ob- 
tain explicit solutions either for the time-dependent case 
with specific initial conditions, nor for a putative steady 
state. However, exact arguments about the behaviour of 
the system may be formulated by mapping the motion 
onto that of an elastic particle. 

We discuss first the semi- infinite case I = oo. The 
following transformation 



r 3 x: 



t^t' = r~H (2) 

leaves the equation of motion (p|) and the variance of the 
noise invariant. Moreover, Eqs (0) imply v — ► v' = r^v, 
so if we perform such a rescaling of variables after each 
collision then in terms of the new variables the motion 
is that of a randomly accelerated particle which collides 
elastically with the wall. We define 



r -3n(t) 



dt 
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where n is the number of collisions. Furthermore, we can 
remove the wall and consider the motion of a free parti- 
cle, in which case n is the number of times the particle 
has crossed the line x = and we should take the ab- 
solute value of x when transforming back to the original 
variables. The bar denotes coordinates referring to the 
fictitious free particle. We can invert Eqn (0) to give: 



t= / dsr 2n ^. 
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We will see that inelastic collapse occurs when n increases 
with t in such a way that t approaches a finite limit. 

To facilitate the discussion, we make a further trans- 
formation onto a stationary Gaussian process (SGP) via 
the following change of variables p3: 
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T = lnt. 



The equation of motion may be written in the form 
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dX 



-X = H(T), 
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where H(T) is a Gaussian white noise with correlator 
(H(Ti)H(T 2 )) = 35(Ti - T 2 ). The correlation func- 
tion for X(T) is (X(T )X(T + T)) = (3/2)e~ T / 2 - 
(l/2)e -3T / 2 . The average number of returns to the ori- 
gin in an interval dT of a SGP with correlation func- 
tion C{T) - 1 - AT 2 at small T is d(n) = pdT, where 
p = V2A/7T Thus d(n) = pdt/t, so (n) = pint for 

large t. In a naive mean-field approach, we would re- 
place n(t) by (n(t)} in Eq. ^ giving t w / 



Therefore, when r < 
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i is Unite. 



This means that after a finite time, the particle has col- 
lided with the wall an infinite number of times and, since 
(x 2 ) = (r 6n x 2 ), where (x 2 ) ~ t 3 , the trajectory collapses 
and the particle adheres to the boundary. 

We now turn to an exact calculation of the return ve- 
locity distribution after many collisions with the bound- 
ary, from which further details of the collapse transition 
can be extracted. In particular, we will show that the 
value of r c predicted by the mean field analysis is correct 
and that the collapse takes place in a finite time due the 
presence of a characteristic velocity scale which decays 
exponentially with the number of bounces. The calcu- 
lation involves two steps. Firstly, we will determine the 
speed distribution on the first return to the origin for 
a particle governed by eqn(Q) and released from x = 
with velocity vq. Secondly, we will use this distribution 
as a Green's function to relate the velocity distribution 
at the boundary after n collisions to that after n — 1 col- 
lisions. By iterating the resulting recursion relation, the 
full return distribution after n bounces can be generated. 

The velocity distribution on the first return to the 
boundary can be calculated from the steady state prop- 
erties of the following escape problem. Particles are in- 
jected into the region x > from the origin at a con- 
stant rate and with initial velocity vq. They experience 
a random force, eqn(^), and can only exit from the re- 
gion at x — 0. The steady state flux of particles leaving 
with speed v is proportional to the first return probabil- 
ity, P(v\vq), for particles released from the origin with 
initial velocity vq. Calculations of this type have been 
carried out in the context of Kramers' equation and 
are known as albedo problems in boundary layer theory. 
Here we will show that the albedo solution for the un- 
damped, random acceleration model with delta-function 
injection has a simple analytic form which, to our knowl- 
edge, has not been noted previously. 

In the steady state, the probability density function 
P{x, v) corresponding to the Langevin equation (Q) obeys 
the Focker-Planck equation 



d 2 P(v,x) dP(v,x) 
= v- 



dv 2 



dx 
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For the escape problem outlined above, this equation 
must be solved subject to the boundary conditions 



P(x, v) — > 0, v — > ±oo 
P(x, v) — > 0, x — > +oo 
P(0,v) = S(v-v ), for v > 



(8) 
(9) 
(10) 



with -P(0, v) for v < the unknown function we wish to 
determine. Using separation of variables, the most gen- 
eral form of the solution to eqn(0) which is continuous 
and differentiable along the line v = and is consistent 
with the boundary conditions eqn((8|) is 



P(x,v) 



d\e- Xx a{\)Ai(-\3v), 
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where Ai is the Airy function and a(A) is as yet unknown. 
Using the orthogonality properties of the Airy functions 
over the the interval [—00,00], this equation can be in- 
verted to express a(A) as an integral over P(x, v). Along 
the line x — this reduces to 

\ia{\) = Ai(-\iv ) + / dvvAi(-\iv)P{0,v), (12) 

J — 00 

where we have used the boundary condition eqn([Io|). The 
requirement that P(x, v) — * for x — * 00 can only be sat- 
isfied if a(A) = for all A < 0. Imposing this condition 
on ( |l2] ) results in an integral equation for the unknown 
part of the boundary velocity distribution, 



Ai(-A3i>o) 



dvvAi(\3v)P(0,-v), 



(13) 



which must hold for all A < 0. 

In the semi-infinite system the only characteristic ve- 
locity is vq, so this necessarily sets the scale of P(0, —v) 
Using the relation P(v\vq) = —vP(0,—v) and the def- 
inition P(v\vq) 
ansatz 



i-/(^), we proceed by making the 
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which was motivated by a numerical simulation of the 
above escape problem. By expressing the Airy function 
in terms of Bessel functions and performing the integral 
in (|l3|), one can verify that P(0,—v) obtained from the 
above ansatz is indeed a solution to ([l3]) for all A < as 
required. 

Next we wish to determine how the return velocity 
distribution evolves as the particle collides many times 
and is reflected with a coefficient of restitution r at each 
bounce. Let P n (v) be the probability density of returning 
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to the wall with speed v after n collisions. This distribu- 
tion obeys the recursion equation 



P, 



W = r (rr) (15) 

n rv' V rv' ) 



n+lV 



with f(x) given by (fL4|) , as particles incident with speed 
v are reflected with speed rv. Eqn(|l5|) may be solved 
by first changing variables to u = bxv which turns 
the integral into a convolution. Using standard Fourier 
techniques, eqn(|l5|), together with the initial condition 
Pq{v) = S(v — i>o), can be solved up to quadrature giving 



Pn(v) 
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gife(ln(u/t) )~nlnr) 



[coshf] 



(16) 



We are now in a position to discuss the collapse tran- 
sition in some detail. Firstly, one can explicitly calculate 
the moments of the velocity distribution P n (v). One finds 



2cosf (1 + a) 



(17) 



for —5/2 < a < 1/2, whilest moments with a outside 
this range do not exist. The fluctuations in v are thus 
extremely large and, as n increases, different moments of 
v diverge for different values of r. This is because v de- 
pends exponentially on n, so the moments are sensitive 
to the rare events associated with the extremes of the 
return velocity distribution. However, if one is interested 
in typical trajectories, the natural variable to consider 
is lnu as this grows only linearly with n fT2[| . Changing 
variables in eqn(fl6|) to u — \tl(v/vq) and defining the nor- 
malized probability distribution Q n {u) = e u P n (e u ), one 
finds that in the large n limit, 



exp 



n In r 



7T 



(18) 



From the peak of the distribution we can identify a crit- 
ical value of r, r c = e _7r /^. Note that this value is 
the same as that predicted by the mean field analysis. 
Furthermore, the distribution of lnu is sharply peaked 
around a value n\n(r/r c ) with fluctuations of order y/n. 
There is thus a typical, characteristic velocity which be- 
haves like (r/r c ) n . On scaling grounds, one would expect 
the time intervals between collisions for r < r c also to de- 
cay exponentially with n since v ~ t 1 / 2 . Consequently, 
the total elapsed time after an infinite number of bounces 
will remain finite and, subsequently, the particle will re- 
main at rest on the boundary. 

We shall now discuss the case of a particle in a fi- 
nite box of size I. The trajectory of a particle starting 
at x — with speed vq is the same as for the semi- 
infinite system up to the instant where x — I for the 
first time, and may be obtained from a trajectory of 
the equivalent SGP and the transformations (||), (||). In 



terms of the variable X(T), the position of the far wall 
is L = Zexp(— |T — 3nlnr). The quantity ln£ decreases 
linearly in T, but increases by — 31nr at each return of 
X to zero. Figure 1 shows the logarithm of the enve- 
lope X nmx of a typical trajectory, defined as the largest 
value of \X\ during the interval between the previous and 
the next zero, obtained by simulating equation (Q), to- 
gether with the corresponding values of lnL(T) for three 
values of r. Because the correlations in the SGP decay 
exponentially in T, the intervals between zeros will have 
short-range correlations. It follows that the motion of 
InL is that of a biased random walker, biased towards 
the origin for r > r c and away from the origin for r < r c . 
Meanwhile, the distribution of \X\ is very sharply cut off 
at \X\ ~ 1. It is known that a random walker with a 
bias away from the origin has a non-zero probability of 
never returning to the origin ||. Thus, for r < r c , the 
particle has a non-zero probability of never reaching the 
far wall. If the particle does reach the far wall, it will 
then have a non-zero probability of never reaching the 
near wall. The particle will therefore collapse to one of 
the walls, and for each trajectory X(T) the time taken 
for the particle to come to rest will be finite. For r > r c 
the particle will always reach the other wall in a finite 
time and the process will repeat indefinitely. 
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FIG. 1. A semilogarithmic plot of the envelope I™, (see 
text) of a typical trajectory of the SGP (Eqn. M), to- 
gether with the effective position of the far wall L(T). When 



r < r c = e 



t/%/3 



0.163..., the particle has a non-zero 



chance of never reaching the wall. 

There is therefore a symmetry-breaking transition be- 
tween a steady state where the particle bounces an infi- 
nite number of times from either wall, and the collapsed 
state where the particle bounces a finite number of times 
from one wall and an infinite number of times from the 
other in a finite time. Figure 2 shows typical trajectories 
for four values of r, obtained from simulations of Eqn (fj) 
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and transforming back to the real variables x and t, show- 
ing the collapse beneath the critical value r c w 0.163. 
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FIG. 2. Typical trajectories of the particle in a finite box 
for four values of r, showing inelastic collapse for r < r c . 



The arguments about collapse of a single particle with 
a wall can trivially be extended to two inelastic parti- 
cles, in which case the collapsed state consists of the 
two particles moving together. Collapse can also be ex- 
pected in a many-body system, when pairs of particles, 
and then larger clusters, aggregate. Real systems may 
deviate from the ideal case studied in this Letter by 
having short-range correlations in the driving force, or 
by the coefficient of restitution approaching unity in the 
small velocity limit. In these cases, one would still expect 
some remnant of the collapse transition, and this should 
provide a method for extracting information about the 
time- and velocity-scales on which such deviations occur. 
We will discuss these and other related points at greater 
length elsewhere [ [l6| . 
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Finally we argue that a collapse transition can occur 
for particles obeying Kramers' equation of motion, 



d 2 x dx , , 



(19) 



This equation is often solved perturbatively by expanding 
about the high viscosity limit [l4|]. However, if collisions 
with a boundary are to be included, the effect of the in- 
ertial term needs to be treated in a more complete way. 
We will again consider a particle colliding with a single 
boundary with coefficient of restitution r. Scaling argu- 
ments analogous to those presented above suggest that a 
collapse transition will take place for all 7 at the same 
value of r as in the undamped case. 

We first perform a rescaling of the variables as in (||). 
The resulting equation of motion has the same form as 
( |l9| ) but with a rescaled 7 given by 7' = r 2 j. Af- 
ter many collisions there is an effective time dependent 
7(F) = r 2n ^7. From Eqn (||), we sees that in the un- 
damped problem for r < r c , r 2n ^ — > faster than l/t. 
If one introduces a time dependent 7 in eqn(|l9|), a simple 
scaling analysis shows that it will be irrelevant asymp- 
totically if it too decays faster than l/t. Thus, the dis- 
sipative system will always be in the collapsed state for 
r < r c . Now let us assume that the dissipative system 
can collapse for some r > r c . As the collapsed state 
is approached, the velocity of the particle goes to zero 
while the acceleration remains of the order of the noise 
strength. Consequently the dissipative term ~/v will be- 
come negligible and, as we have assumed r > r c , the 
particle will ultimately move away from the wall. We 
thus conclude that the collapse transition will occur when 
r = r c = g-^/v's f or anv value of 7 fl5|]. 
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